Thermal analysis of production of resonances in relativistic heavy-ion collisions 
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Production of resonances is considered in the framework of the single-freeze-out model of ultra- 
relativistic heavy ion collisions. The formalism involves the virial expansion, where the probability 
to form a resonance in a two-body channel is proportional to the derivative of the phase-shift with 
respect to the invariant mass. The thermal model incorporates longitudinal and transverse flow, 
as well as kinematic cuts of the STAR experiment at RHIC. We find that the shape of the 7r + 7r~ 
spectral line qualitatively reproduces the preliminary experimental data when the position of the 
p peak is lowered. This confirms the need to include the medium effects in the description of the 
RHIC data. We also analyze the transverse-momentum spectra of p, K*(892)°, and /o(980), and 
find that the slopes agree with the observed values. Predictions are made for r\, vj ', u>, <j>, A(1520), 
and E(1385). 
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I. INTRODUCTION 



With the help of the mixed-event technique the STAR 
Collaboration has managed to obtain preliminary results 
on the production of hadronic resonances. The data were 
presented for IT (892)° p(770)° H0, / (980) SB, 

and A(1520) HI10- 

The invariant-mass distribution of 
the particles produced in the decays of the resonances 
may be used to get the information on the possible in- 
medium modifications of hadron masses. Such modifi- 
cations are predicted by various theoretical models of 
dense and hot hadronic matter 0, |?j , as well as hinted by 
the enhancement of the dilepton production in the low- 
mass region [HI [D| • Most interestingly, the measured 
invariant-mass distribution of the 7r + 7r~ pairs 0, 0] also 
suggests a drop of the effective mass of the p meson by 
several tens of MeV 3j. This effect was recently dis- 
cussed by Shuryak and Brown 0, Kolb and Prakash 
0, and Rapp 0, with the conclusion drawn that it 
is a genuine dynamical effect induced by the interaction 
of the p meson with the hadronic matter. On the other 
hand, the measurement of the invariant- mass distribu- 
tions of the ttK pairs indicates that the effective mass of 
the if* (892)° remains unchanged 0- 

Since one measures the properties of stable hadrons 
which move freely to detectors, the experimentally ob- 
tained correlations may bring us information only about 
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the final stages of the evolution of the hadronic mat- 
ter, i.e., about the conditions at freeze-out. From this 
point of view it is interesting to analyze the production 
of resonances in the framework of the thermal approach 
(the single-freeze-out model of Refs. [l5l jig, Il7j ) which 
turned out to be very successful in reproducing the yields 
and spectra of stable hadrons. In particular, the thermal 
model can be naturally used to study the impact of the 
possible in-medium modification of the p-meson mass on 
different physical observables, e.g., on the correlation in 
the invariant mass of the 7r + 7r - pairs, the ratios of the res- 
onance abundances, or their transverse- momentum spec- 
tra. 

The paper is organized as follows: In the next Sec- 
tion we outline the Dashen-Ma-Bernstein formalism used 
to describe a gas of hadronic resonances. In Sec. 11111 
we present the invariant-mass correlations of 7r + 7r - pairs 
emitted from a static thermal source, and in Sec. II VI we 
discuss the effect of the temperature of such a source 
on the shape of the spectral 7r + 7r - line. In Sec. [3 
we present the main assumptions of the single-freeze-out 
model. With the knowledge of both the experimental 
kinematic cuts (Sec. IVI|) and the feeding from the decays 
of higher resonances (Sec. IVII|) . the model is then used 
to compute the invariant mass spectrum of the 7r + 7r - 
pairs (Sec. IVIII|) . In Sec. IIXI and IXl we present the 
model results for the ratios of the resonance yields and 
for the resonance transverse-momentum spectra (wher- 
ever it is possible we compare the model results with the 
preliminary data). The paper contains three Appendices 
which give the parameters for the 7r + 7r~ phase shifts and 
explain in simple terms the implementation of the exper- 
imental kinematic cuts in our calculations of two- and 
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three-body decays. 



II. PRODUCTION OF RESONANCES AND 
THE PHASE-SHIFTS 

The formalism for the treatment of resonances in a 
hadronic gas in thermal equilibrium has been developed 
by Dashen, Ma, and Bernstein [T^ . and Dashen and Ra- 
jaraman 19], and in the context of heavy- ion physics 
has been recalled and further elaborated by Weinhold, 
Friman, and Norenberg jU HJ [H |H ( see also 
[25l Efil 127^ ) . The basic formula following from the for- 
malism is that density of the resonance per volume and 
per unit invariant mass, M, produced in the two-body 
channel of particles 1 and 2 in thermal equilibrium is 
given by the formula 



dn 
dM 



f 



d 3 p 1 dS 12 (M) 



1 



(2n) 3 



dM 



exp 



T 



(1) 



± 1 



where S12 (M) is the phase shift for the scattering of par- 
ticles 1 and 2, / is a spin-isospin factor, T is the tempera- 
ture, and the sign in the distribution function depends on 
the statistics. In practice, one may replace the distribu- 
tion function by the Boltzmann factor, since the effects 
of the quantum statistics are numerically small. 

As pointed out by the authors of Ref. 22], in many 
works the spectral function of the resonance is used ad 
hoc as the weight in Eq. (QJ instead of the derivative of 
the phase shift, which is the correct thing to do. For 
narrow resonances this does not make a difference, since 
then both the spectral function and the derivative of the 
phase shift peak very sharply at the resonance position, 
ra#, i.e., dSi2(M)/dM ~ ttS(M — ra#), and then one 
recovers the narrow-resonance limit 
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(27T) 



exp 



(2) 



± 1 



However, for wide resonances, or for effects of tails, the 
difference between the correct formula Q and the one 
with the spectral function is very significant, not only 
conceptually but also numerically. 

We first focus on the 7r + 7r~ scattering. We use the 
experimental phase shifts, which can be parameterized 
with simple formulas of Ref. |23| : 



tan (5f(M)) = y^l 



ml/M 2 q 21 
Ctf + D{<f){Aml 



(3) 



M 2 



where q = \J M 2 - 4mf/2, = 139.57 MeV is the 
mass of the charged pion, and the remaining parameters 
are listed in App. E] The relevant channels are 1=1, 
I = 1 0), I = 0, I = (fo/a), and I = 2,1 = 0. 



Figure [I] shows the weights in Eq. JI}, i.e., the quan- 
tities fd8{M)/{ndM), plotted as functions of M. The 
factor / consists of the spin degeneracy, 21 + 1, and the 
isospin Clebsch-Gordan coefficient, equal 2/3 for the isos- 
inglet, 1 for the isovector, and 1/2 for the 1 = 2 channel. 
In the isovector channel the peak of the weight function 
is found at 762.7 MeV, which corresponds to the inflec- 
tion point in S(M), and not the point where S(M) = 7r/2, 
which is at 773.6 MeV (note a 10 MeV "dropping" of the 
p- meson "mass"). For the isoscalar channel the corre- 
sponding curve in Fig. Speaks at 932.5 MeV. We note 
that in this channel the weight function fd5/(irdM) is 
significantly different from the spectral function. The lat- 
ter one appears in the calculation of the production rate 
(see for instance 0, whereas our study concerns 

the invariant- mass spectra, where formula (1) holds. As 
a result of the use of the phase-shifts, there is a very small 
contribution from the a meson, since the strength of the 
scalar channel, as seen in Fig. ^ is small. It is worth- 
while to notice that it peaks at the two-pion threshold 
[2^ . which is an immediate consequence of Eq. © for 
1 = 0. The contribution of the 1 = 2 channel is tiny 
(also note that it is negative since the phase shift in this 
channel is a decreasing function of M, cf. PJj). 

For three-body decays formulas analogous to Eq. JQ 
may be given, cf. Ref. HIDil. They would involve the 
detailed dynamical information on the decay process. For 
our purpose this is not necessary (the three-body decays 
feed the range of rather low invariant 7r + 7r~ masses) and 
also not practical, since the detailed information on the 
dynamical dependences of the appropriate transition ma- 
trices is not easy to extract from the experimental data. 
In our calculations we include the resonances decaying 
into three-body channels which are very narrow (a;, 77, 
and 7/). Thus, as is customary in similar applications, 
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FIG. 1: The weight functions, 1/irf d5/dM ', in the isoscalar, 
isovector, and isotensor channels, plotted as a function of the 
7r + 7r - invariant mass. The phase shifts are taken from the 
parameterization of the experimental data |28|. Numerical 
labels indicate the peak positions (in MeV) which correspond 
to inflection points in the phase shifts as functions of M. The 
1 = 2 channel is negligible. 
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we only account for the phase-space dependence on the 
invariant mass M, and do not consider dynamical effects 
of the transition matrix. 



III. RESONANCES IN THE STATIC SOURCE 

First, in order to gain some experience, we consider the 
simplest case of the emission of particles from a static 
source. We will come back to a more realistic descrip- 
tion incorporating the flow in Sec.|V| We also assume in 
this Section the single- freeze-out hypothesis 0, 0, ^3 ? 
hence no effects of rescattering are incorporated after the 
chemical freeze-out. We thus assume the chemical freeze- 
out temperature to be [2^, l3l| 



^ chem 



165 MeV, 



(4) 



and, as we have said, the thermal freeze-out occurs at the 
same temperature, T t herm = ^chem- In the next Section 
we will lift this assumption. 

We now simply use Eq. with T = T c h e m, and in- 
clude the following resonances that couple to the 7r + 7r - 
channel: p, fo/cr, rj, r/ , Ks, and /2(1270). For the 
omega, both the three-body mode, u — > 7r + 7r _ 7r°, and the 
two-body mode, u — > 7r + 7r _ , are incorporated. The ap- 
propriate branching ratios are included in the constant / 
in Eq. QJ. For channels other than those of Fig.d where 
the experimental phase shifts are included, we use the 
simple Breit-Wigner parameterization, which is good for 
the narrow resonances. For the width of Ks we take the 
experimental resolution of the STAR experiment, which 
is about 10 MeV [23. The width of the oj is also increased 
by the same value. 

We compute the spectra at mid-rapidity, hence we use 



dn 



dMdy 



y=0 



i J p l r 



p±dp± 1 dSi(M) 
(2tt) 2 n dM 



(5) 



■Pi 



exp 



where i labels the channel, n is the volume density of the 
7T + 7r - pairs, and M is the invariant mass of the 7r + 7r~ 
system. The above formula, written for two-body decays, 
is supplemented in the actual calculation with the three- 
body reactions. The limits for the transverse- momentum 
integration are taken at the upper and lower cuts for 
the STAR experiment 4], p 1 ™ = 0.2 GeV and p h ^ gh = 
2.2 GeV. 

The results of the calculation are shown in Fig. 12(a). 
The two-body resonances lead to well- visible peaks, while 
the three-body channels of the lj, r\ and 7/ produce typical 
broad structures at lower values of M. Note that the 
scalar-isoscalar channel generates the resonance, /o, and 
a smooth background, a. This background increases with 
the decreasing M, and peaks at the threshold. 



IV. CHEMICAL VS. KINETIC FREEZE-OUT 

In this section we analyze in simple terms the effects of 
possible lower temperature of the thermal freeze-out on 
the shape of the 7r + 7r - invariant-mass correlations. To 
this end we simply use Eq. © at a lower value of T. 

Consider, for brevity of notation, the gas consisting 
only of 7r + , 7T~, and p° . We assume isospin symmetry for 
the system, which works very well for RHIC The 
pions rescatter elastically through the p-channel. The 
total number of pions that are observed in the detectors 
is equal to 



N„ = r hem (2x^(T chem ,/i = 0) 



+ 2 x 3 x n p (T chem ,/i = 0)) ; 



where 



d 3 p 
(2tt) 3 



■ exp 



-i/mf +p 2 + ii 



(6) 



(7) 



The factor of 2 in front of in Eq. © comes from the 
7r + -7r - degeneracy, the factor of 2 in front of n p comes 
from the fact that rho decays into two pions, and the 
factor of 3 is the spin degeneracy of the p. The chemical 
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FIG. 2: The volume density of the n tt~ pairs plotted as 
a function of the invariant mass, obtained from the thermal 
model with a static source, Eq. (0. (a) The case of Ttherm = 
T c hem = 165 MeV. (b) The case of T the rm = HO MeV. The 
numbers indicate the position of the p peak (in MeV). The 
labels indicate the decays included, with the two-body tt + tt~ 
channels and the three-body channels, to — > 7r°7r + 7r _ , 77 — > 
7r°7r + 7r - , and 7/ — > r]7T + 7T~. The relative contributions from 
the resonances are fixed by the value of the temperature. 
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potential, /i, vanishes at the chemical freeze-out. The 
volume at the chemical freeze-out is denoted by y chem . 

The chemical freeze-out is defined as the stage in the 
evolution when the abundances of stable (with respect to 
the strong interaction) particles have been fixed. Next, 
the system expands, and while it progresses in its evo- 
lution, the elastic scattering may occur. This scattering 
proceeds through the formation of unstable resonances, 
until the system is too dilute for the scattering to be ef- 
fective and the thermal freeze-out occurs. At the moment 
of the thermal freeze-out the total number of pions in our 
example is 

JV,r = ^ therm (2xn^(T therm ,^) 

+ 2 x 3 x n p (T therm , fip)) , (8) 

where y therm is the volume at the thermal freeze-out, 
and lip = (the reaction p° <-> 7r + 7r~ is in equilib- 
rium) are the chemical potentials which ensure that N n 
is conserved, as requested by the fact that the system 
had frozen chemically. Knowing the ratio y therm /y chem 
we could compute p n comparing the right-hand sides of 
Eqs. © and ®. However, for the present purpose, where 
we are only interested in the shape of the correlation 
function in M and not the absolute values, this is not 
necessary. The chemical potential enters in Eq. ((HJ as 
a multiplicative constant, V therm exp(2/i 7r /T t herm)- Since 
we do not control in the present calculation the volume, 
we can remove the normalization constant from our con- 
sideration. The argument holds if more reaction channels 
are present. 

Thus, we redo the calculation of the previous Section, 
based on Eq. ((5J), but now with a lower temperature 
and an arbitrary normalization constant. The results 
for Ttherm = 110 MeV are shown in Fig. (b). We 
note a prominent difference from the case of Fig. [3| (a), 
where the temperature was significantly higher: the high- 
M spectrum is suppressed, while the low-M spectrum is 
enhanced. Beginning from the high-mass end, we notice 
that the ji resonance has practically disappeared, the 
relative strength of the fo to the p peak is only about 
1/5 compared to 1/3 in Fig. 0(a), finally, left to the K s 
peak we note a significantly larger background from the 
a tail, the shoulder of the p, and the r] decays. While 
for T = 165 MeV the hight of this background is only 
slightly above the hight of the p peak, for T = 110 MeV 
it rises to about two times the hight of the p peak. 

We also remark that due to the presence of the thermal 
function in Eq. (0, the position of the peak is shifted 
downwards from the original vacuum value to 752.2 MeV 
for Ttherm = 165 MeV, and to 744.0 MeV for T th erm = 
110 MeV. 

To conclude this Section, we state that the very sim- 
ple thermal analysis shows that the shape of the "spec- 
tral line" of the 7r + 7r~ system depends strongly on the 
temperature of the thermal freeze-out, which, when com- 
pared to accurate data, may be used to determine T t herm 
in an independent manner. In the next sections we elabo- 



rate on this observation by incorporating other important 
effects. 



V. FLOW AND THE SINGLE-FREEZE-OUT 
MODEL 

The medium produced at mid-rapidity in ultra- 
relativistic heavy-ion collisions undergoes a rapid expan- 
sion, characterized by the longitudinal and transverse 
flow. Although flow has no effect on the invariant mass 
of a pair of particles produced in a resonance decay, since 
the quantity is Lor entz- invariant, it nevertheless affects 
the results, since the kinematic cuts imposed in the ex- 
periment in an obvious manner break this invar iance. In 
Ref. 0, we have constructed a thermal model which 
includes the flow effects. The model has the following 
main ingredients: 

1. There is one universal freeze-out, occurring at the 
temperature 

T chem = T therm = T ~ 165 MeV. (9) 

in El El El III we have demonstrated that with 
the inclusion of all hadronic resonances the distinc- 
tion between the traditionally considered chemical 
and thermal freeze-outs is not necessary, at least 
for RHIC. 

2. To describe the geometry and flow at the freeze-out 
we adopt the approach of Refs. [51 |3 El El El 
El El liol IZlL |42 | . The freeze-out hypersurface is 
defined by the condition 



r = yjt 2 — r\ — r 2 — r 2 z = const, (10) 

while the transverse size of the fire-cylinder is made 
finite by requesting that 

P = \/ r l+ r l <Pmax. (11) 

In addition, we assume that the four- velocity of the 
expansion at freeze-out is proportional to the coor- 
dinate, namely 

^ t ( r x r y r z \ . N 

uM = V = 7( 1 'f'f'f)- (12) 

The model is explicitly boost-invariant, which in 
view of the recent data delivered by BRAHMS E! 
is justified for the description of particle production 
in the rapidity range — 1 < y < 1. 

3. All resonances from the Particle Data Tables E3 
are included in the calculation, which is important 
due to the Hagedorn-like behavior of the resonance 
mass spectrum H! El El E3 



5 



The model has altogether four parameters: two ther- 
mal parameters, T, and the baryon chemical potential, 
/ie, which are fitted to the available particle ratios, and 
two geometric parameters, r and p m ax, fitted to the 
transverse-momentum spectra. The details of the model 
can be found in Ref. 17]. The model works very well 
and economically for the particle ratios [2lj, transverse- 
momentum spectra [T^ as we ll as strange-particle pro- 
duction [l6j . 



VI. KINEMATIC CUTS 

The next important effect is related to the experimen- 
tal cuts and traditionally has been the domain of ex- 
perimentalists. However, in the present application the 
inclusion of all relevant kinematic cuts of the STAR anal- 
ysis 0, m> LiJ can be included straightforwardly. Needless 
to say that the proper inclusion of kinematic constraints 
is frequently crucial when comparing theoretical models 
to the data. 

For two-body decays, the relevant formula for the num- 
ber of pairs of particles 1 and 2, derived in App. [0 nas 
the form 



dN 12 _ d5 12 bm f p ^sh ± 



/•Pi.high m,hi g h rPhigh rvhigh 
/ dp 1 / dy x dp dy 

Jptlow ^ 2/1, low Jyiov, 



dM dM p* 

I sin 70 1 

with all quantities defined in App. |Bl 

For the case of three-body decays we have 



(13) 



12 



dM 



POO PC 

) / 27rp- L dp 1 - 
Jo J - 



dN 12 (p ± ,y) ± 
d V 7T7 S (P )' 



dM 



(14) 



where b is the branching ratio and dNi 2 (p ± ,y)/dM has 
been evaluated in App. with full inclusion of all cuts 
relevant in the STAR experiment. 

The cuts in the STAR analysis of the^^Tr^ invariant- 
mass spectra have the following form 



W < i, 

1^1 < 0.8, 
p^ > 0.2 GeV, 



(15) 



while the bins in pt = + Pn I start from the range 
0.2 - 0.4 GeV, and step up by 0.2 GeV until the range 
2.2 — 2.4 GeV. These conditions are implemented in the 
calculation shown below. 



VII. DECAYS OF HIGHER RESONANCES 

An important ingredient and, in fact, the key to the 
success of the thermal models in both reproducing the 



particle ratios and the transverse-momentum spectra 
is the inclusion of resonances. Although 
the thermal distribution suppresses the high-mass parti- 
cles, their abundance grows exponentially according to 
the Hagedorn hypothesis [45|, |46|, H3, |4g , and in practice 
one needs to go very high up in the mass of the reso- 
nances in order to obtain stable results We include 
all resonances from the Particle Data Tables |44]. The 
high-lying resonances decay in cascades, eventually pro- 
ducing stable particles. 

The resonances considered in this paper, in particu- 
lar the p, also acquire substantial contributions from the 
higher resonances, e.g., rj' — > p°7, or 0(1020) — > pir. Such 
effects, entering at the level of a few tens of %, are diffi- 
cult to account for accurately in our formalism. This is 
due to the dynamics characterized by the parameters that 
are not well known. For instance, the decay of rj f may 
proceed through the p 0/ y channel, as well as directly into 
the uncorrelated three-pion state, 7r + 7r _ 7r°, which forms 
a smooth background around the p peak. In other words, 
the feeding of the p from the higher resonance need not 
reproduce the shape of the p peak from Fig. 1, and the re- 
sulting dependence on M may be altered to some extent. 
This important but, due to experimental uncertainties, 
not easy issue is left for later studies. At the moment we 
take the simplest way, and assume that the shape of the 
spectral line in Fig. 1 is not altered. This amounts to 
including a multiplicative factor, d, for each considered 
resonance. These factors are obtained from the thermal 
model as discussed in Ref. [17], |29( . The thermal param- 
eters for the full RHIC energy of ^snn = 200 GeV are 



T 

PB 
PS 
Pi 



165.6 MeV, 
28.5 MeV, 
6.9 MeV, 
-0.9 MeV. 



(16) 



The calculation leads to the following enhancement fac- 
tors coming from the decays of higher resonances: dx s = 



1.98, d n 



1.74, d a 



1.13, d p = 1.42, d u 



1.43, 



drj> = 1.08, d fo = 1.01, and d h = 1.28. Thus, the effects 
is strongest for light particles, Ks, rj, p, and lj, while it 
is weaker for the heavier rj f and scalar mesons. 



VIII. RESULTS OF THE SINGLE-FREEZE-OUT 
MODEL 

We may finally come to presenting the results of 
the calculation in the framework of the single-freeze-out 
model. The calculation includes the kinematic cuts de- 
scribed in App. El an d O an d the enhancement factors 
from the higher resonances, described in Sec. IVIII 
The expansion parameters are taken to be r = 5 fm and 
Pmax = 4.2 fm, which according to the fits of the p± spec- 
tra, corresponds to the centrality 40-80% [HJ. In Fig. [31 
we show the results obtained with the help of Eqs. f|13!14|) 
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FIG. 3: The invariant 7r + 7r - mass spectra in the single- freeze-out model for four sample bins in the transverse momentum of 
the pair, pt, plotted as a function of M. The contribution of various resonances are labeled in Fig. (c), with n indicating the 
joint contribution of n and rj' ' . The kinematic cuts of the STAR experiment are incorporated. The labels at the p peak indicate 
its position in MeV, lowered as compared to the vacuum value due to thermal effects. 



with the STAR kinematic cuts (|15|). for four sample bins 
in the transverse momentum of the pair, p^, with the 
lowest pt in Fig. 3(a) and highest in Fig. 3(d). The con- 
tributions of various resonances are clearly visible. We 
note that the shape of the spectrum changes with the as- 
sumed pt bin. For both the lowest and the highest pr the 
low-M region is suppressed, while at intermediate pr the 
spectrum has large contributions at low M. In our calcu- 
lation the relative contributions form various resonances 
is fixed, as given by the model. Also, the relative nor- 
malization between Figs, (a-d) is preserved. The labels 
at the p peak indicate its position in MeV. We observe 
that it is lowered compared to the vacuum value, due to 
the thermal effects, and assumes values between 740 and 
747 MeV. This fall is not as large as that observed in the 
preliminary STAR data 0, Q , where the position of the 
peak resulting from fitting the data is much lower, about 
700 MeV. Our single-freeze-out model with the vacuum 
p meson is not capable of producing this result. 

This is perhaps the most important outcome of our 
analysis: both the naive thermal approach, cf. Fig. 2, 
and the full-fledged single-freeze-out model with expan- 
sion and kinematic cuts, are unable to reproduce the (pre- 
liminary) STAR data if the vacuum value of the 1=1 



7T7T phase shift (3J) is used. This confirms the earlier con- 
clusions of Shuryak and Brown and Rapp ^4(. The 
model with the vacuum p can bring the peak down to 
~ 740 MeV, calling for about additional 40 MeV from 
other effects, such as the medium modifications. 

In order to show how the medium modifications will 
show up in the 7r + 7r~ spectrum, we have scaled the tttt 
phase shift in the p channel, according to the simple law 

^i(M) scaled = ^ 1 (5- 1 M) vacuum , (17) 

where s is the scale factor. We use s = 0.91, which places 
the peak at 700 MeV. The result of this calculation is 
shown in Fig. Q] Now, with the lowered p position, the 
calculation looks very similar to the preliminary data of 
Ref. 0, • With the p peak moved to the left the dip 
between the p and Ks peaks is largely reduced, as seen 
in the preliminary data. Also, the background on the left 
side of the Ks peak, coming from other resonances, is, 
for the considered pr bin, significantly higher than the 
hight of the p peak. 

We also remark that the shape of the the curves ob- 
tained in the single-freeze-out model for intermediate 
values of pr is closer to the naive thermal-model cal- 
culation of Fig. 2 with Ttherm = 110 MeV rather than 
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FIG. 4: Same as Fig. Etc) with the p meson scaled down by 
9%. Now the peak of the p has moved to the preliminary 
value observed in the STAR experiment 



the masses of other resonances (except for the mass of 
if* (892)° which is not changed). In Ref. gg we studied 
the effect of the in-medium mass and width modifications 
on the outcome of the thermal analysis in the situation 
where a common scaling of baryon and meson masses 
with the temperature or the density was assumed (only 
the masses of the pseudo-Goldstone bosons were kept 
constant). The results of Ref. 49] showed that moderate 
modification particle masses are admissible. A satisfac- 
tory description of the ratios of hadron abundances mea- 
sured at CERN SPS may be obtained in a thermal ap- 
proach with the modified masses of hadronic resonances, 
however, the changes of the masses affect the values of 
the optimum thermodynamic parameters. Similar results 
were obtained from the analysis of the first RHIC data 
The effects of the in-medium mass and width mod- 
ifications on the outcome of the thermal analysis were 
also studied in Refs. [H El El El El 



^therm = 165 MeV. This is reminiscent of the "cooling" 
effect of the resonance decays, cf Ref. 29]. Indeed, 
the enhancement factors di, higher for low mass M and 
lower for high mass M, connected with the feeding of 
the resonances by even higher excited states, lead to the 
modification of the spectral line such that it resembles a 
lower temperature spectrum obtained without the reso- 
nance feeding. 



IX. RATIOS OF THE RESONANCE YIELDS 

In this and the next Section we present the results of 
the single-freeze-out model for the ratios of the resonance 
yields and the resonance transverse-momentum spectra. 
In our analysis we shall consider two cases: in the first 
(standard) case we assume that the masses of hadrons 
are not changed by the in- medium effects, whereas in the 
second case we assume that the mass of the p meson at 
freeze-out is smaller than its vacuum mass. Inspired by 
the preliminary STAR results H Q , we shall use for this 
purpose a rather low value of m* = 700 MeV. In these 
two cases the values of the thermodynamic parameters 
are determined from the experimental ratios of the yields 
of stable hadrons (those listed in Table 1 of Ref. [5l|). 
This gives T = 165.6±4.5 MeV and p B = 28.5±3.7 MeV 
for the standard case |5l|, and T = 167.6 ±4.6 MeV and 
Pb — 28.9 dz 3.8 MeV for the case with the modified mass 
of the p meson. The reason for not including the ratios 
of the resonance yields as the input in our appraoch is 
twofold: firstly, the data describing the production of the 
resonances are preliminary; secondly, there are no mea- 
surements of the resonance spectra in the central events. 
For example, the data on p° production were collected 
for the centrality class 40-80%. 

We take into consideration the modification of the 
mass of the p meson only, since at the moment there 
are no experimental hints concerning the behaviour of 



Our predicitions for the ratios of the resonance yields 
are presented in Table 1. The most interesting are the 
results for the p°/tt~ ratio which is 0.11 for the case with- 
out the in- medium modifications and 0.14 for the case 
with the modified rho mass. The first value is in a good 
agreement with a theoretical number presented recently 
by Rapp [l4]. We can see that the effect of the dropping 
p-meson mass helps us to get closer to the preliminary ex- 
perimental result of 0.18, however the theoretical and the 
experimental values still differ by more than one standard 
deviation. Our model value for the ratio fo (982) /tt~ is 
a factor of four smaller than central value of the prelimi- 
nary experimental result. On the other hand, our results 
for K*(892)/K~ and A(1520)/A are larger than the pre- 
liminary experimental values, whereas the *(892) ra- 
tio agrees rather well with the experiment (cf Table 1). 

A comparison of the theoretical and experimental re- 
sults for the ratios involving resonances is interesting 
from the point of view of the discussion on the decou- 
pling temperature, T t herm, ^.e., the temperature char- 
acterizing the thermal freeze-out. If T t he . r ™ is signifi- 
cantly smaller than T c h em one expects (53, EH E3 that 
the measured ratios such as K* / K or p°/tt~ are smaller 
than the values determined at the chemical freeze-out. 
This behavior is due to the readjustment of the reso- 
nance abundances (formed in elastic collisions) to the 
decreasing temperature on the path the system follows 
from T c hem down to T t herm- Our results shown in Ta- 
ble 1 indicate that certain experimental ratios are in- 
deed smaller (such as K*(S92)/K~ or A(1520)/A for 
central collisions), whereas some are larger (p°/ir~ and 
/o(982)/7T~). Consequently, from the comparison of the 
theoretical and experimental ratios, at the moment, it is 
hard to draw a definite conclusion about the difference 
between the two freeze-outs. The study gets even more 
involved in view of possible in-medium modifications of 
the masses, which affect the ratios and may complicate 
the overall thermodynamic picture. 
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Model 
m* = 770 MeV 


Model 
m * = 700 MeV 


Experiment 


T [MeV] 


T= 165.6 ±4.5 


T= 167.6 ±4.6 




p B [MeV] 


Ms = 28.5 ±3.7 


p B = 28.9 ±3.8 




77/71-" 


0.120 ±0.001 


0.112 ±0.001 






0.114 ±0.002 


0.135 ±0.001 


0.183 ±0.028 [4] (40-80%) 


(jj/7T~ 


0.108 ±0.002 


0.102 ±0.002 




i^*(892)/7r- 


0.057 ±0.002 


0.054 ± 0.002 






0.025 ±0.001 


0.024 ± 0.001 






0.0121 ±0.0004 


0.0115 ±0.0003 




/o(980)/tt- 


0.0102 ±0.0003 


0.0097 ± 0.0003 


0.042 ± 0.021 [4] (40-80%) 


iT(892)/if- 


0.33 ±0.01 


0.33 ±0.01 


0.205 ±0.033 [4] (0-10%) 
0.219 ±0.040 [4] (10-30%) 
0.255 ± 0.046 [4] (30-50%) 
0.269 ± 0.047 [4] (50-80%) 


<I>/K* (892) 


0.446 ± 0.003 


0.448 ± 0.002 


0.595 ±0.123 [4] (0-10%) 
0.633 ±0.138 [4] (10-30%) 
0.584 ±0.132 [4] (30-50%) 
0.528 ±0.106 [4] (50-80%) 


A(1520)/A 


0.061 ±0.002 


0.062 ± 0.002 


0.022 ± 0.010 [7] (0-7%) 
0.025 ± 0.021 [7] (40-60%) 
0.062 ± 0.027 [7] (60-80%) 


E(1385)/E 


0.484 ± 0.004 


0.485 ± 0.004 





TABLE I: Predictions of the thermal model for the ratios of hadronic resonances measured at RHIC at ^/snn — 200 GeV. 
In the two model calculations, the two different effective masses of the p meson have been assumed. In the first case (second 
column) m* = 770 MeV, whereas in the second case (third column) m* = 700 MeV. The two thermodynamic parameters 
have been fitted only to the ratios of stable hadrons. In the first case they are the same as those discussed in Sec. 2. The 
preliminary experimental data are taken from 4, 7]. Both the experimental and the theoretical pion yields are corrected for 
the weak decays. 



X. TRANSVERSE-MOMENTUM SPECTRA 

In this Section we present the single-freeze-out model 
results for the transverse-momentum spectra of various 
hadronic resonances. The method of the calculation is 
the same as that used in the calculation of the spectra of 
stable hadrons 0,^30- I n particular, the feeding of 
the resonance states from all known higher excited states 
is included, which leads to the enhancement of the res- 
onance production characterized by the factors d{. As 
usual, the two thermodynamic parameters and the two 
geometric parameters (fitted separately for different cen- 
trality windows) are taken from Ref. [5l|. Knowing the 
centrality dependence of r and p ma x we may analyze the 
resonance production at different centralities and com- 
pare it to the existing data (note that the data on p pro- 
duction are collected only for rather peripheral collisions 
corresponding to the centrality window 40-80%). 

In Fig. we show our results for the transverse- 
momentum spectra of p° and /o(892), and compare them 
to the preliminary data obtained by the STAR Collabo- 
ration Q. The expansion parameters are r = 5 fm and 
Pmax =4.2 fm, which corresponds to the centrality of 40- 



80% 51]. In the presented calculation the vacuum value 
of the 7T7T phase shift has been used. We have checked 
that the scaled phase shift gives very similar results, and 
the small change of the thermodynamic parameters due 
to the drop of the rho mass affects the spectrum very lit- 
tle. The results shown in Fig. [5| indicate that our model 
can quite well reproduce the experimental spectrum of 
the rho meson and the slope of the /o(980) spectrum. We 
find, however, the discrepancy in the normalization be- 
tween the model and the preliminary experimental data 
for the /o(980) spectrum, which reflects the result of Ta- 
ble 1. In Fig. El we show our predictions for the spec- 
tra of A(1520) and £(1350). Since A(1520) is currently 
measured in the central collisions 7], in this case we 
used the values of the geometric parameters correspond- 
ing to the most central collisions, namely r = 8.7 fm 
and Pmax = 7.2 fm, which corresponds to centrality of 
— 20% 51]. In Fig. |H|we also show the spectrum of 
£(1385) which might be measured by STAR in the near 
future p]. In Fig. 0we compare our predictions for the 
spectrum of if * (892) to the preliminary experimental re- 
sults [l|. To complete our discussion, we also show our 
predictions for other resonances for the most-central case, 
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FIG. 5: The mid-rapidity transverse- momentum spectra of 
the p° and f (980) mesons, as obtained from the single-freeze- 
out model. The model parameters, r = 5 fm and p ma x = 
4.2 fm, correspond to centralities 40 — 80 % 51]. The data 
points are taken from Ref. The vacuum value of the tttt 
phase shift is used in the model calculation, with the scaled 
phase shift yielding very similar results. 
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FIG. 6: The single- freeze-out model predictions for the mid- 
rapidity transverse-mass spectra of A(1520) and £(1385). The 
model parameters, r = 8.7 fm and p ma x = 7.2 fm, correspond 
to the most central events, c = — 20% [HJ. 



namely 77, p, cj, </>, and 7/ and /q. 



XI. CONCLUSION 

In this paper we have presented a comprehensive anal- 
ysis of resonance production at RHIC in the thermal ap- 
proach. We have studied the 7r + 7r - invariant mass spec- 
tra, the ratios of the resonances, and their transverse- 
momentum spectra. 
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FIG. 7: The mid-rapidity transverse- momentum spectra of 
the various resonances, as obtained from the single- freeze- 
out model. The preliminary data for the K* are taken from 
Ref. [l|. The model parameters, r — 8.7 fm and p ma x — 
7.2 fm, correspond to centralities — 20 % |51|. 



Our calculation of the invariant mass distribution of 
the 7T + 7r - pairs, performed in the framework of the single- 
freeze-out model, has been done with the use of the 
derivatives of the experimental phase shifts (and not the 
spectral function), which guarantees the thermodynamic 
consistency of our approach. Decays of higher resonances 
have been included in a complete fashion. Moreover, 
we have taken carefully into account the realistic kine- 
matic cuts corresponding to the STAR experimental ac- 
ceptance. We find that the preliminary STAR result in- 
dicating a drop of the rho meson mass down to 700 MeV 
cannot be reproduced in our model when the vacuum 
properties of the p are used. A good qualitative agree- 
ment is found, however, when the mass of the p is low- 
ered. Therefore our calculation confirms the dropping- 
mass scenario for the p, induced by the medium effects. 

We stress that the shape of the 7r + 7r~ spectrum is very 
sensitive to the freeze-out temperature. In a sense, it 
can be used as a "termo meter" , independent of other 
methods of determining the temperature, such as the 
study of the ratios of the stable particles, or the slopes of 
the transverse momenta. We have shown that at small 
T the contributions from the low lying resonances are 
enhanced, and the contributions at high mass are sup- 
pressed. On the contrary, at large T the contributions 
from the high lying states are enhanced and become visi- 
ble (like, e.g., the contribution from /2(1270)), while the 
low mass contributions are reduced. In the single-freeze- 
out model, although the freeze-out process takes place 
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at relatively high temperature T ~ 165 MeV, the de- 
cays of the resonances lead to an effective "cooling" of 
the spectrum, with low-mass resonances acquiring more 
feeding from the higher states than the high-mass res- 
onances. This would be equivalent to taking a lower 
temperature in a naive calculation without the resonance 
decays. Such effects, seen initially in the shape of the 
transverse-momentum spectra, are therefore also seen in 
the invariant mass 7r + 7r - distributions. 

The experimental information on the resonance pro- 
duction is crucial for our understanding of the freeze-out 
process. In particular, it may be used to distinguish be- 
tween the scenario with the two freeze-out s, separate for 
the elastic and inelastic processes, and the scenario with 
a single freeze-out. For example, in the scenario with the 
two distinct freeze-outs the ratios involving the resonace 
yield in the numerator and the stable-hadron yield in the 
denominator (e.g., p/n, K*/K, or A(1520)/A) should be 
smaller than the analogous ratio calculated at the chemi- 
cal freeze-out. By comparing the preliminary STAR data 
to the predictions of the thermal model we have pointed 
out that certain ratios of this type are indeed somewhat 
smaller, while other are larger. Consequently, at the 
present time one cannot exclude which approximation, 
the one with two or the one with one freeze-out, is more 
appropriate in the thermal approach. More acurate data 
would be highly desired in that regard. 
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APPENDIX A: PARAMETERS FOR THE tttt 
PHASE SHIFTS 

The following parameterization is used for the tttt 
phase shifts of Eq. (J3J 12a] (all in units of = 
139.57 MeV): 



4° 
^0 



0.220, A\ = 0.379-10" 
0.268, B\ = 0.140- 10" 4 , B { 



A 2 ■ 

2 




-0.444 • 10" 
-0.857-10 



r 2 



-0.0139, C\ = -0.673 • 10 
-0.00139, £>1 = 0.163- 10" 
36.77, s{ = 30.72, s 2 = -21.62. 



7 D 2 



-0.00221, 
-0.129- 10~ 3 , 
(Al) 



APPENDIX B: KINEMATIC CUTS FOR 
TWO-BODY DECAYS 

Let the unlabeled quantities refer to the decaying res- 
onance, and labels 1 and 2 to the decay products. The 
number of pairs coming from the decays of the resonance 
formed on the freeze-out hypersurface is 0, 



TV- 



/ lt Cl J ^CC 2 B( P , Pl )S(p), (Bl) 

where the symbols C, C\ , and C2 denote kinematic cuts, 
and the source function of the decaying; resonance is ob- 
tained from the Cooper- Frye formula ED] > 



S(p) = J (x) p" f\p-u(x)]. 



(B2) 



The quantity / is the thermal distribution of the particle 
decaying at the hypersurface E, with the collective flow 
described by the four- velocity u. We pass to rapidity and 
transverse momentum variables in the laboratory frame, 
and have explicitly 

"<>■») <B3) 

b „ / m^mj; cosh(?/ — yi) 



p ± Pi cos 7 



m 



E* 



where pi and E* denote the momentum and energy of 
particle 1 in the rest frame of the decaying resonance, 
7 is the angle between p 1 - and p± in the laboratory 
frame, and b is the branching ratio for the considered 
decay channel. 

In general, for the considered cylindrical symmetry, we 
have 



C 



61 blow < p < Phi g h)%iow < y < 2/high), 



(B4) 



d = °(pt\ow < Pi < ^hi g h) 6> (^,iow < y < 2/i,high), 

with i = 1, 2 and 9 (a < x < b) — 1 if the condition is sat- 
isfied, and otherwise. Due to momentum conservation 
we should use in (|B4|) 



2/2 



= sinh -1 



+ (pi) — 2p ± p^- cos 7 
m 1 - sinh?/ — sinh?/i 



mi 



(B5) 



Next, we perform the integration over 7 in Eq. (|B1IB3|) . 
The 5 function gives the condition cos 7 = cos 70 , with 

m ± mi cosh(y - y ± ) - mE{ 
cos 70 = — r — l , (B6) 



P^Pi 



which takes effect only if — 1 < cos 70 < 1. A factor of 
2 follows from the two solutions of 7 = arccos(cos7o) for 
7 e [0, 2tt). The final result is 



1777, /*Pl,high 
T / d 



dpi 



yi, high 



brn 

Pl J ^r,iow J yi, low 
, 0(-l < cos 70 < 1) 



dyi 



sin 70 



S(p^ 



Phigh rVhi S h 
dp / dy 



(B7) 
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where we have introduced 



C 1 ® C 1 I 

^2 2lcos7=cos7o i 



(B8) 



with the substitution (|B5|) understood. 

The experimental cuts frequently involve cuts on 
pseudo-rapidity of particles. This amount to adding ex- 
tra conditions of the form 



??- ow <^lo, 




m^ 2 sinh 2 ^ + ra^sinh^ 



ra^sinh 2 ^ — mj- 2 smh. z y 



(B9) 



which may be implemented in Eq. (|B7|) in the form of 
functions, denoted by C] 7 . 

The above formulas have been written for a sharp res- 
onance. For a wide resonance, according to Eq. Q, the 
generalization assumes the form (JT3|) . 



APPENDIX C: KINEMATIC CUTS FOR 
THREE-BODY DECAYS 

The three-body decay is usually considered in the rest 
frame of the resonance of mass m. Here, since the cuts 
are defined in the laboratory frame, we need to consider 
the kinematics in this frame. The phase-space integral 
for the decay of particle of mass m and momentum p 
into products 1,2, and 3 is proportional to 



/ 



d 3 pi d 3 p 2 d 3 p 3 
Ei E2 E3 
X(5 (3) (p-pi -p 2 



C5{E-E 1 -E 2 -E 3 ) 
-Ps)\M\\ 



(CI) 



where we have introduced E = yj m 2 + p 2 , e£c, and C 
denotes a generic cut in the kinematic variables, to be 
specified later. For simplicity we assume that Ai can be 
approximated by a constant, i.e. only the phase-space ef- 
fect is included. This condition can be relaxed at no dif- 
ficulty, if needed. We are interested in the invariant- mass 
distribution of particles 1 and 2, hence we introduce the 
S[M - ^(E 1 + E 2 ) 2 - (pi +p 2 ) 2 ] function in Eq. (IHT1) 
and obtain the expression for the probability of emitting 
a pair of invariant mass M from a particle moving with 
momentum p: 



dN l2 (p) 
dM 



d 3 pi d 3 p 2 d 3 p 3 



CS(E-E 1 -E 2 -E 3 ) 

5 (M - V(£i + £ 2 ) 2 -(Pi+P 2 ) 2 ) 

5^ (p - pi - p 2 - P3 ) , (C2) 



where a is a normalization constant. First, we trivially 
carry the integration over p 3 through the use of the last S 
in Eq. (|C2|) . Next, we pass to the rapidity and transverse- 
momentum variables, and carry the integration over the 



angle between momenta p^ and p^, denoted as 7. This 
yields 



dN 12 (p^,y) 
dM 



= 27ra^ J Pidp^dyida J p 2 dp2dy 2 C 
1 M 

x -A^^T- \ S (A-E 3 ), (C3) 

Apfpf I sin 70 1 



where a is the angle between p 1 - and pi, 

-M 2 + m 2 + m\ + 2m^m2"Cosh(yi — y 2 ) 



cos 70 



(C4) 



and the sum over e results from the two solutions for 70, 
i.e. sin 70 = eyl + cos 2 70. Next, 

A = m^coshy — mi coshyi — m^cosh^, (C5) 
and, finally, 

E 3 = [B-2p ± pi cos a (C6) 

1/2 

-2p ± p 2 L cos(a - 70) + 2p^p2 cos 70] 

with 



B — m 3 + (m sinh?/ — m 1 sinh?/i — rn 2 sinh^) 2 



+p ±2 +pi 2 +pi 2 - 



(C7) 



We need still to carry the integration over the angle 
a. We square the expression under the S function in 
Eq. (|C3|) . obtaining 

A 2 = B-2p ± pi cos a (C8) 
— 2p ± p 2 L (cos a cos 70 — sin a sin 70) + 2pipi cos 70. 



The squaring imposes the condition 

A > 0. 



(C9) 



This equation can be solved straightforwardly by intro- 
ducing t = tan(a/2), 



cos a 



sm a = 



1-t 2 
1 + t 2 ' 
2t 

TTt 2 ' 



(CIO) 



and the notation 



G 



-A 2 -\- B -\- 2 cos J0P1P2 



P^P2 



rr . Pi 
H COS 70 H r . 

P2 

Now Eq. (|C6|) acquires the simple quadratic form 



(Cll) 



(1 + t 2 )G = 2te a/1 - cos 2 70 + (1 - t 2 )H, (C12) 
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with the solutions 



to(e^e) = tan(ao(e, e')/2) 



(C13) 



with e' 
dition 



_ ey/l - cos 2 70 + e' \/ H 2 - G 2 + sin 2 70 
G + H ' 

±1. The solutions make sense under the con- 



H 2 -G 2 + sin 2 70 > 0. 



(C14) 



From the derivative of the delta function we obtain the 
factor 



A 



P^-P2 \H sinao(e, e') — e sin 70 cosao(e, e')\ 



(C15) 



It is easy to check that this factor is independent of e and 
e' . Thus, to the extent that the cut function C does not 
involve azimuthal angles, we may use one combination of 
these signs and put a factor of four. The final result is 



dN 12 (p ± ,y) 
dM 



Sira / dp 



?idyi J dp 2 



dy 2 C 



x 0(A)0(H 2 -G 2 + sin 7 2 ) 
M 1 



(C16) 



I sin 70 1 p^P2 I H sin ao — sin 70 cos ao | 



where ao is any of the angles (|C13|h and all the nec- 
essary substitutions are understood. The normalization 
constant a can be obtained from the condition 



/ 



dM dN "%/ V) = 1, 
dM 



(C17) 



with no cuts present, 
unity, (7=1. 



i.e. with the cut function set to 



The form of the cut function C involves the ranges in 
the integration variables p± , y\ , p^ , 2/2 ? the cuts on the 
pseudo-rapidity of particles 1 and 2, as well as cuts on 
the rapidity and the transverse momentum of the pair of 
particles 1 and 2 [23 . These cuts assume a simple form of 
products of the 6 functions. Then, Monte Carlo methods 
are appropriate to compute Eq. (|C16|) . 
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